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We study the accuracy of the pion form factor, obtained with a local-duality version of dispersive 
sum rules. To probe this accuracy, we make use of a potential model, where the exact form 
factor may be calculated from the solution of the Schrodinger equation and confronted with the 
local-duality form factor. The deviation between these quantities is found to be below 20% in 
the region of momentum transfers Q > 2 — 3 GeV, independently of the specific form of the 
confining potential. We argue that the local-duality model for elastic form factors in QCD has at 
least this level of accuracy. 
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1. Introduction 

The pion elastic form factor at Q = 3 — 8 GeV is sensitive to the specific features of the onset of 
the perturbative regime and opens the possibility to study subtle details of the pion structure. The 
experimental study of the pion form factor at the lower part of this region will become available 
with an upgrade of JLab in the next few years. Therefore it is now the right time to review the 
theoretical understanding of this quantity. Several versions of the method of QCD sum rules may 
and have been applied to this problem. However, no conclusive results have been obtained and we 
still have a strong discrepancy between the results from different versions of QCD sum rules [|TJ, 
The goal of the presented work was to study the accuracy of the pion form factor obtained from 
so-called local-duality (LD) version of QCD sum rules Let us briefly remind the basic ideas 
leading to the LD model. 



2. Sum rule 

The basic objects for sum-rule calculations of pion properties are the correlation functions 

n(p 2 ) = f (a\Tj{x)f(o)\a)^dx, 
r{pi P 2 2 ,q 2 )= J{a\Tj(x 1 )j(o)j\x 2 )\a)e i P x ^dx 1 dx 2 . (2.1) 

Here Q. is the physical vacuum; j(x) is a short-hand notation for the interpolating axial current 
jsa(x) °f the positively charged pion, (Q. [j5a(0)| 7i(p)) = ip a f n ', J{0) denotes the electromag- 
netic current J v (0). For brevity, we omit Lorentz indices. In QCD these correlators may be cal- 
culated by applying the operator product expansion (OPE). Instead of the Green functions in the 



Minkowski space ( |2.1| ), it is convenient the evolution operators in the Euclidean space, which 
emerge after performing the Borel transform p 2 — > X, the parameter z being related to the Eu- 
clidean time. The Borel transform leads to several improvements: (i) suppresses the contributions 
of the excited states; (ii) improves the convergence of the perturbative expansion; (iii) provides the 
necessary smearing required by quark-hadron duality. The Borel image of the two-point correlator 
has the form 

f'OO 

rioPE^) = J p pe n(s)e- st ds + U cond (z), ppert(^) = Po(s) + a s pi{s) + 0(a 2 ). (2.2) 

Here p,-(s) are the spectral densities of the two-point diagrams of the perturbation theory, Yl CO nd(t) 
describes nonperturbative power corrections. Making use of the hadron intermediate states, for the 
two-point correlator we obtain 



U(z) = f 2 e m * T + excited states. (2.3) 

The first term in this expression corresponds to the pion contribution. 

The double Borel transform p\ 2 — > | of the three-point function has the form 

r°° r°° _ s j+fi 
r O PE(T,0= / / A per tOi,j 2 , Q)e ^ t ds 1 ds 2 + r cond (i:,Q), 
Jo Jo 

Apert(si) S2, Q) = Aq{sus 2 , Q) + c&AiCsi, s 2 , Q) + 0{a 2 ), (2.4) 
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Apert(si , S2, Q) is the double spectral density of the three-point diagrams of the perturbation theory. 
r c ond(^i Q) describes power corrections. Making use of the hadron intermediate states gives 

r(T, Q) = F n {Q)fle~ m * % + excited states. (2.5) 

The key assumption of the method of sum rules is the duality assumption which says that the 
contribution of the excited states is dual to the high-energy region of the perturbative diagrams. 
Using this assumption, the sum rules take the form in the chiral limit of the vanishing quark masses 



(2.6) 



si = fV(^ + ^ + ^l^ + ---, 

Jo 12 % 81 

Fn{Q)fl = / / A pen (s i ,s 2 ,Q)e ^ z ds x ds 2 

Jo Jo 

(f G 2 ) 4%a s (qq) 2 , , , , 

+ 24 T + ll^ T ( 13 + G T ) +- " (2 " 7) 

These relations are the standard relations for the extraction of the hadron decay constants and form 
factors in the method of QCD sum rules. 



Let us focus on Eq. (2.7). We would like to study the form factor at large Q. The form factor of 



a bound state should decrease with Q; however, the power collections of the r.h.s. are polynomials 



in Q and thus rise with Q. So, Eq. (2.7) cannot be directly used at large Q. There are two ways 
for considering the region of large Q: The first way is the resummation of power corrections, after 
which the resummed power correction decrease with Q. This may be done in a model-dependent 
way by making use of nonlocal condensates [0J. The second way is just to set the Borel parameter 
T = 0; then all power corrections vanish and the remaining perturbative part decreases with Q. This 
version of sum rules is called a local-duality (LD) sum rule [||]. In the LD limit one finds 

fl = rPvMds = ^ (l + ^) + O (a, 2 ) , (2.8) 



o 



4 k 2 V K 

f*eff(2) /-ieff(e) 



F*{Q)fl = J j o A pen (s l ,s 2 ,Q)ds l ds 2 . (2.9) 

The double spectral densities p P ert(s) and A pert (5 1 , s 2 , Q) are given by the perturbation theory; f n is 
known from the experiments. So, if we fix s e ff(<2), the form factor may be calculated. The spectral 
densities have the following properties: at Q — > the spectral densities of two- and three-point 
functions are related to each other by the Ward identity 

IimAi(ji,^,fi) = p i (si)S{s 1 -s 2 ). (2.10) 
At Q — > oo explicit calculations give: 



1 8 K 

lim AoOi,j 2 , Q) ~ T^T, lim Ai(j b j 2 , Q) = t^Po(>i)Po(>2)- (2.11) 



For the pion form factor, two rigorous properties are known: The normalization condition related 
to the current conservation 

F„(0) = 1. (2.12) 
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The asymptotic behavior at large Q due to the QCD factorization theorem 

Q 



F n {Q) = ^4^ + -- (2-13) 



Obviously, if we set 

s^(0 -> 0) = ^ 



s eS (Q 0) = , s eS (Q -> oo) = 4n 2 f 2 , (2.14) 
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then the form factor obtained from the LD sum rule ( |2.9| ) satisfies both of these rigorous properties. 



The two values of the s e s at large and at small Q are not far from each other! So, it is easy to 



construct an interpolation function s e ff(2) f° r a U Q with the limiting values (2.14). 
Now, we can formulate the LD model for hadron elastic form factors: 1 

(i) It is based on a dispersive three-point sum rule at % = (i.e. infinitely large Borel mass 
parameter). In this case all power corrections vanish and the details of the non-perturbative 
dynamics are hidden in a single quantity - the effective threshold s e ff(<2)- 

(ii) It makes use of a model for s e ff (2) based on a smooth interpolation between its values at 
Q — > determined by the Ward identity and at Q — > oo determined by factorization. Since 
these values are not far from each other, one believes the details of this interpolation to be 
not essential. For instance, a self-consistent expression may be used [JXJ] : 

4tt 2 f 2 

'-08) = — ^ (2-15) 

Thus, the only non-perturbative input for the LD model is the pion decay constant f n . 

Obviously, the LD model is an approximate model which does not take into account the details 
of the confinement dynamics, and it is important to understand its accuracy. Now, where this 
accuracy may be tested? 

The only property of theory relevant for this model is factorization of hard form factors. There- 
fore, the model may be tested in quantum mechanics for the case of the potential containing both 
the Coulomb and the confining interactions. The spectral representation for the form factor and the 
decay constant in the LD limit are similar to those in QCD; the corresponding spectral densities 
can be calculated from the two- and three-point diagrams of the non-relativistic field theory. To 
probe the sensitivity of the LD model to the details of the confining potential, we shall consider 
two different confining potentials 

(X m (O 2 r 2 

V{r) = --+V conf {r), V conf (r) = or (1) V conf {r) = — ^— (2) (2.16) 

and make use of the parameters relevant for hadron physics. The parameters are chosen such that 
the Schrodinger equation for both confining potentials leads to the same value of decay constant 
(i.e. »F(r = 0)): m = 0.35 GeV, co = 0.5 GeV, a = 0.168 GeV, a = 0.3. The exact form 
factors as obtained from the solution of the Schrodinger equation are different in these two models; 
however the LD model for the form factor, which depends only on the value of a and *F(r = 0) for 
both models is the same. Comparing the exact form factors and the LD form factors allows us to 
probe the accuracy of the LD model. 

'Notice, however, that the model is not expected to work at small nonzero Q as the OPE is not applicable here. 
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3. Numerical results 

Fig. [j] presents our results for the potential model. The black lines give the results from the 
LD model for the effective threshold and the corresponding form factor. The blue and the red lines 
are the exact form factors obtained from the solution of the Schrodinger equation (Fig. Plot: la) and 



the corresponding exact thresholds which reproduce these form factors by the LD expression (2.5). 
[The variable k is related to the variable s, used above as s = 4(k 2 + m 2 )]. Recall that the form 
factor in our potential model behaves as F(Q 2 ) ~ l/Q 4 because we consider spinless quarks. The 
lesson to be learnt from the potential model is the following: the exact threshold fc e ff does not exceed 
the LD threshold by more than 5% As Q increases, the accuracy of the LD approximation increases 
rather fast, too. This conclusion does not depend on the details of the confining interaction. 

Fig. ^| shows the analysis of the pion form factor in QCD. The black lines - the results from 
the LD model; the red lines represent our upper boundary for the threshold and the corresponding 
form factor. According to the experience from the potential model we set the upper boundary for 
the threshold * e ff by about 10% higher than the asymptotic threshold. We believe this estimate 
to be quite reliable: the previous analysis of various correlators [||, ^] shows that the extraction 
procedures in quantum mechanics and in QCD are very similar both qualitatively and quantita- 
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Figure 1: (a) The exact vs LD form factors, (b) The exact vs LD effective thresholds. Red lines - linear 
confining potential; blue lines - harmonic oscillator confining potential; black lines - LD model. 
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Figure 2: (a) Models for the effective threshold, (b) The corresponding pion form factor. Black lines - 
the LD model. Red lines - our upper boundary. Green lines - the results from the sum rule with non-local 
condensates [g]. Data from [gfl. 
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tively. Notice also that the LD model for the pion form factor agrees well with the results from the 
dispersion approach 

The problem is that the effective threshold recalculated in |Q] to reproduce the results from the 
sum rule with non-local condensates exceeds our upper boundary by more than 10%. This might 
be a consequence of adopting in Ref . [Q] a procedure of fixing a T-independent effective threshold 
based on the maximal stability. According to the results of [^] this algorithm does not guarantee 
the extraction of reliable values. 

4. Summary and Conclusions 

We studied the LD model for the elastic form factor which may be formulated in any theory 
where the form factor at large momentum transfers satisfies the factorization theorem (i.e., any 
theory containing both Coloumb and Confining interactions). 

Our main conclusions are: 

1. In the region 2=1—2 GeV, the exact effective threshold exhibits rapid variation with Q. 
Depending on the theory (relativistic or non-relativistic), the error of the LD form factor in this 
region may reach 30-40% level. In general, for a relativistic theory a smaller error is expected. 

2. At Q > 2 — 3 GeV, the LD model provides a good description of the pion elastic form factor - 
with the accuracy better than 20%. Moreover, the accuracy increases rather fast with Q. 

We point out that our prediction for the pion form factor is considerably lower than the pre- 
diction of the approach based on the sum rule with non-local condensates. This discrepancy needs 
clarification. Presumably [^], its origin might be traced back to the procedure of fixing the z- 
independent effective threshold in the method of sum rules with non-local condensates based on 
merely the Borel stabilty criterion. We are going to redo the analysis making use of the recently 
formulated modifications of the sum-rule method based on the T-dependent effective threshold [0]. 
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